Probing Quantum Phase Transitions on a Spin Chain with a Double Quantum Dot 
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Quantum phase transitions (QPTs) in qubit systems are known to produce singularities in the 
entanglement, which could in turn be used to probe the QPT. Current proposals to measure the 
entanglement are challenging however, because of their nonlocal nature. Here we show that a double 
quantum dot coupled locally to a spin chain provides an alternative and efficient probe of QPTs. We 
propose an experiment to observe a QPT in a triple dot, based on the well-known singlet projection 
technique. 
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Spin chains have been studied for many years because 
of their simple formulation, which enables analytical so- 
lutions, and their similarity to more complex quantum 
many-body systems. Recently, it has been possible to en- 
gineer fully tunable spin chains based on quantum dots 
with one or more electrons pQ. Such chains can be used 
as spin qubits for quantum computing [2], or as a spin 
bus whose ground state transmits quantum information 
over large distances [5HS]- In all these settings, the ability 
to mediate entanglement is paramount for incorporating 
spin chains into quantum devices. 

Quantum phase transitions (QPTs) can have a strong 
effect on adiabatic operations involving the ground state 
of a spin chain. QPTs occur at energy level crossings as 
a function of external parameters, between ground states 
with very different physical properties [7j . In the critical 
regime where the two ground states are nearly degener- 
ate, macroscopic observables such as two-qubit entangle- 
ment can exhibit non-analytic behavior [8-12J. In finite- 
size systems, the phase transition is typically discontin- 
uous, or first order; however, the underlying physics is 
essentially the same as in infinite systems. 

From a quantum information perspective, QPTs may 
produce singularities in the entanglement, which could 
potentially enhance device operation. Alternately, we 
could view entanglement as a sensitive probe of the 
ground state, which could potentially enhance our un- 
derstanding of QPTs. The latter approach was recently 
applied to XY-type spin chains, using nonlocal pairs of 
qubits to probe the QPT [T5HT5] . In this configuration, 
strong enhancements of the entanglement were observed 
near critical points. However, in realistic quantum dots, 
the spin couplings are not of the XY type, and nonlocal 
measurements can be rather challenging [16 . 

Because of its simplicity, a local probe could poten- 
tially be more effective. Unfortunately, the simplest type 
of measurement - mapping out the effective interaction 
between the spin chain and a weakly coupled qubit - does 
not exhibit unusual behavior near an energy level cross- 
ing [17] . It has been suggested that the time evolution of 
a coupled qubit could be used to probe a QPT [TSJ [T5] 



via the Loshmidt echo [20]. This method has been suc- 
cessfully implemented using nuclear magnetic resonance 
(NMR) [2TJ [22] . Such dynamical methods are powerful, 
but they are still more challenging than simple projective 
measurements. 

In this Letter, we propose a local, projective scheme 
for detecting QPTs in a quantum dot spin chain. The 
chain undergoes successive QPTs as a function of the ex- 
ternal magnetic field. We study two different external 
probes of the QPT. First, we consider a nonlocal probe, 
consisting of two qubits weakly coupled to the chain at 
different locations. We calculate the entanglement be- 
tween the probe qubits numerically, using the concur- 
rence measure [23], and we observe singularities when 
the spin chain undergoes a QPT. We also obtain analyt- 
ical estimates for the entanglement using perturbation 
theory. Next, we consider a local probe, consisting of a 
double quantum dot coupled to a single node of the spin 
chain. We find that the ground state properties of the 
chain are imprinted onto the probe, and we investigate 
the concurrence singularities both numerically and ana- 
lytically. Interestingly, we find that the probability for 
the probe qubits to form a singlet state echos the non- 
analytic response of the concurrence. This is significant 
because the singlet probability is relatively easy to mea- 
sure, using the singlet projection technique common to 
spin qubit experiments |24j . We propose a simple exper- 
iment to test these concepts on the smallest possible spin 
system of size N — 1. This "chain" has a single energy 
level crossing as a function of magnetic field, and the 
ground state transition exhibits a non-analyticity consis- 
tent with a QPT. Our proposal involves a total of three 
quantum dots, and it is therefore within reach of current 
triple dot technologies [25H27] , 

Energy level crossings in the spin chain. — We adopt 
an isotropic Heisenberg model of a spin chain, as ap- 
propriate for single-electron spins in quantum dots. The 
Hamiltonian for a chain of length N is given by H c = 
Jc^fj^i s j ' s j+i — B c ^2j = x Sj z where Sj are spin op- 
erators for the individual electrons. The bare exchange 
couplings between the spins are labelled J c , and the ap- 
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FIG. 1. Graphical representations of the Hamiltonians stud- 
ied here, (a)-(c) correspond to a nonlocal probe in which 
the probe qubits (1 and 2) are attached to different nodes on 
the spin chain, (d)-(f) correspond to a local probe in which 
a pair of qubits is attached to a single node on the chain. 

(a) and (d) describe the full, physical geometry, where the 
chain (lightly shaded circles) is formed of an arbitrary num- 
ber of physical spins. Here, we show probe qubits attached 
to the endpoints of the chain; however, similar results are 
obtained for any attachment points, (b), (c), (e), and (f) 
represent effective geometries, in which the spin chain in its 
ground state is replaced by a pseudospin, as indicated by a 
filled circle (when appropriate), (b) and (e) represent non- 
critical Hamiltonians (far away from a critical point), where 
the bus ground state is non-degenerate (pseudospin-O). (c) 
and (f) represent critical Hamiltonians, where the bus ground 
state is doubly-degenerate (pseudospin-1/2). In the effective 
Hamiltonians, the bus pseudospin interacts with the probe 
qubits via effective couplings (J) and effective fields (73). In 

(b) , the effective coupling J12 is weak (i.e., second order), as 
indicated by a dashed line. 
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FIG. 2. (a), (b) Energy spectra of spin chains of length TV = 4 
and 5, with no coupled qubits, as a function of the magnetic 
field B c expressed in energy units, and scaled by the coupling 
constant J c . The energy level crossings of the ground state 
are indicated by circles, (c), (d) The corresponding concur- 
rence between the probe qubits, when they are coupled in the 
geometry of Fig. [l^a). Here, C is dimensionless, and we take 
Jic=J2c~0.02J c . The singular features of the concurrence 
occur at energy level crossings of the chain. 



plied magnetic field is B c = B c i. Throughout this paper, 
we will adopt J c as the unit of energy. Magnetic fields 
will also be expressed in energy units. A pictorial repre- 
sentation of the chain Hamiltonian is given by the lightly 
shaded circles in Figs. [ija) and (d). For now, we ignore 
any couplings to external qubits, and calculate the en- 
ergy spectrum for H c as a function of B c . The results 
are shown in Figs. [2]ja) and (b) for the cases N — 4 and 
5. In this letter, we are most interested in the two lowest 
energy levels for a given B c , whose crossings are indicated 
by circles in Fig. [5] Each level crossing is associated with 
a QPT in the finite-size spin chain. 

Nonlocal probe. — Next, we consider external qubits, la- 
belled 1 and 2 in Figs, [lja) and (d), which will serve as 
probes of the QPT. The probe coupling Hamiltonian is 
given by H p = Ji c Si • si + J2c&2 ■ sjv where are the 
spin operators for the probe qubits. The probes may be 
coupled to any node of the spin chain, with similar re- 
sults. For definiteness here, we have attached them to 
the endpoints of the chain. 

When the probe couplings are turned on, the energy 
levels of the chain expand into energy manifolds. Our 
goal is to probe the QPT without disturbing it, so the 
manifold structure in Figs, ^[a) and (b) should remain 



largely undisturbed. This places constraints on the probe 
couplings. First, the bare coupling constants must be 
small, such that J\ c , J2 C *C J c - The magnetic field ap- 
plied to the probes should also be much smaller than 
J c , necessitating a magnetic field gradient between the 
qubits and the chain. For definiteness, we take the mag- 
netic field on the probe qubits to be zero. Although large 
field gradients are difficult to achieve in the laboratory, 
we will focus on QPTs occurring at zero field, in an odd- 
size chain. For this case, the field gradient is small, and 
does not pose a serious experimental challenge. 

Figures [2jc) and (d) show the concurrence between 
the probe qubits. The calculation is performed after first 
tracing out the spin-chain degrees of freedom from the 
full Hamiltonian, H = H c + H p , to obtain the reduced, 
bipartite density matrix. The concurrence exhibits singu- 
larities which are correlated with the energy level cross- 
ings of the spin chain, as expected for QPTs. Away from 
the level crossings, the concurrence falls quickly to zero. 
The only exception to this regular behavior is observed 
near zero field for even-size chains, where the concurrence 
plateaus at its maximum value, C = 1. 

This interesting behavior can be understood intuitively 



3 



by treating the probe qubits as a perturbation. It has 
previously been shown that when the bare qubit-chain 
coupling is small, and when the chain is in its ground 
state, the interactions can be described by an effective 
Hamiltonian [T7]- Far away from a QPT, the system is 
noncritical and the effective Hamiltonian involves only 
the external spins, as shown in Fig. [ljb): 



N=4 



N=5 



H nc — BiSi z + B 2 S 2z 



E 

a—x,y,z 



Jl2aSi a S2 a 



(1) 



Here, the bare coupling parameters J\ c and J 2c are hid- 
den inside_the effective coupling J 12 and the effective lo- 
cal fields -Bi,2- Note that the effective couplings are gen- 
erally anisotropic, except in special cases. If the qubit 
couplings are turned on adiabatically, the chain will re- 
main in an inert, effective pseudospin-0 state. The ef- 
fective coupling arises due to virtual excitations of the 
chain outside its ground-state manifold; it is therefore 
second order in the perturbation: J\ 2 /J c ~ {Jic/Jc) 2 - 
On the other hand, the effective fields emerge at first 
order: B i2 /J c ~ (Jic,2c/Jc) ■ (Analytical expressions 
for both quantities are provided in |T7j and [55].) We 
therefore generally find that B\_ 2 J\ 2 in the noncrit- 
ical regime. Hence, the external qubits align with the 
effective field to form a separable state, for which C ~ 0. 
The only exception is the special case near B c = for 
an even-size chain. Here B 12 =0 due to the spin-singlet 
character of the chain ground state. Since J\ 2 ^ 0, it 
can generate maximal entanglement between the probe 
qubits, as indicated in Fig.[2^c). 

The situation is very different near a QPT. In this case, 
the ground state of the chain is approximately two-fold 
degenerate and behaves as an effective pseudospin-1/2, 
as shown in Fig. [ljc). The effective Hamiltonian at the 
critical point then describes a simple three-body system: 



lea S\ a S ca + J2ca.S2a.Sca. 

a=x,y,z 



(2) 

Here, the spin operator S c acts on the the pseudospin-1/2 
of the chain ground state. In contrast with the noncrit- 
ical regime, the effective couplings are now first order in 
the perturbation: Ji Ct2c / J c ~ J\c,2c/ Jc- Because J\ c _ 2c 
are relatively large, H cp can mediate entanglement be- 
tween the two probe qubits, with the resulting value of 
the concurrence determined by the relative size of J\ c _ 2c 
compared to B\ 2 . The couplings J\ c .2c enhance entan- 
glement while the fields B\ 2 suppress it. 

We can take this analysis further for the QPT occur- 
ring at B c — in an odd-size spin chain [e.g., the central 
peak in Fig{2](d)]. In this case, the effective Hamiltonian 
has a much simpler, isotropic^form, with Ji c =J 2c =J, 
01,2=0, B C =B C , and H cp = J(Si + S 2 ) • S c — B c S CiZ . 
We can derive an expression for the concurrence between 
the probe qubits which shows the explicit form of the 
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FIG. 3. The concurrence C and singlet probability Ps (dashed 
black and solid red curves, respectively) of a double-dot, local 
probe coupled to a spin chain of size (a) TV = 4 and (b) N — 5, 
as the chain undergoes QPTs. Here we take Ji2=J2 C =0.02 J c . 
C and Ps are both dimensionless. 



non-analyticity [28] : 



C = 



1 - 



\B n 



J/2 



Bl + lJ* + J\B c 



(3) 



We see that the concurrence attains a maximum value 
of 1/3 at the critical point [25], and decreases to zero as 
(J c /B c ) 2 when B c > J c . The full-width-at-half-max of 
the peak is given by (-1 + v/32/5)J ~ 1.53J. 

We close this section by describing an experimental 
procedure for observing concurrence peaks associated 
with QPTs in a spin chain. The spin system is prepared 
in its ground state via thermalization, or an adiabatic ini- 
tialization procedure. The concurrence measurement re- 
quires performing full quantum state tomography of the 
two-qubit reduced density matrix. This involves 15 sepa- 
rate measurements of the two-qubit correlators {cr^o^j}, 
where an (o 2 j) is a Pauli operator acting on qubit 1 (2), 
with i, j e {I, X, Y, Z} [3D]. (We exclude the trivial two- 
qubit identity operator.) 

Local probe. — The local probe geometry that we con- 
sider is shown in Fig. [lid). Here, one side of a double- 
quantum dot is attached to one node of a spin chain. 
In this case, the probe Hamiltonian is given by H p = 
J12S1 -S2 + J 2c S 2 - Si. Using the methods described above, 
we can compute the concurrence between the spins in the 
probe double dot, obtaining the results shown in Fig. [3] 
for two different size chains. We find that the QPTs oc- 
curring in the chain are imprinted onto the probe. In 
Fig. [3] we also show the overlap probability Ps between 
the probe qubits and a singlet state. We see that Ps 
mirrors the singularity in C . 

We can understand the main features in Fig. [3] by ap- 
plying perturbation theory to the local probe geometry. 
As before, the chain is effectively inert away from a crit- 
ical point, and the noncritical effective Hamiltonian is 
given by 



B 2 S 2 , Z + J12S1 • S 2 



(4) 
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FIG. 4. Solid red curves show the singlet probabilities Ps for 
a double quantum dot coupled locally to a spin chain of length 
N, as a function of the magnetic field on the chain. From top 
to bottom, the curves correspond to N — 9,7,5,3, and 1. 
The dashed black curve shows the corresponding concurrence 
C for the case N = 1. Here, we take Ji2=J2 C =0.02J c . 



The entanglement between the probe qubits is deter- 
mined by the interplay between J12 which enhances the 
concurrence, and B2 which suppresses it. The effective 
local field, B2 = (0|si z |0), depends only on the true spin 
of the ground state of the spin chain, |0). In contrast with 
the pseudospin, which remains fixed at 0, the true spin 
increases by 1 in each successive noncritical region, away 
from B c = 0. Accordingly, the concurrence is suppressed 
in discrete steps. 

Approaching a critical point, the chain becomes 
pseudospin- 1/2, and the effective Hamiltonian becomes 

Hep — H nc — B C S CZ 4- S J2caS2aS ca . (5) 

a=x,y ,z 

The singularities all have a downward-pointing "valley" 
shape, which can be explained by considering the B c =0 
transition in Fig.^b). In the non-critical region, Bi ^ 0, 
while in the critical region, B2—O (the latter is only 
true for the B c =0 transition [21]). This would nor- 
mally lead to an upward-pointing singularity, since B2 
suppresses the entanglement. However, a second effec- 
tive coupling (J2 C ) emerges in the critical region, as in- 
dicated in Eq. |5]), which reduces the entanglement be- 
tween qubits 1 and 2 by sharing the entanglement with 
pseudospin c. We now show that the latter effect is al- 
ways dominant, leading to valley-type singularities, by 
obtaining exact solutions for the asymptotic critical and 
noncritical behaviors. 

We focus strictly on the B c =0 transitions in odd- 
size spin chains. Figure [4] shows several B c =0 transi- 
tions for chains of varying length. At the special point 
B C =Q, we find that B c =B2—0, and the effective cou- 
pling J2C is isotropic, yielding the effective Hamiltonian 
^cp = (J12S1 4- </2cS c ) • S2. Since C and Pg are dimen- 
sionless quantities, at the point B c =0 they can only be 
functions of the dimensionless ratio 7 = J2cl J\i- From 
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FIG. 5. (a) Depth of C and Ps valley singularities, like those 
shown in Fig. [4] (with N = 1), when two probe qubits (1 and 
2) are coupled to a quantum dot (c), plotted as a function of 
the coupling ratio. The inset shows the confinement profile 
of the triple-dot experiment, (b) Full-width-at-half-minimum 
(FWHM) of the valleys with J 12 fixed at 0.02J C . In both 
panels, the concurrence C and the singlet probability Ps are 
represented by dashed black and solid red curves, respectively. 

Eq. Q we obtain [25] 

c =l( 1+ 7i=h)- ft -d +<«>/«■<•> 

The concurrence C was also obtained numerically for this 
special case in Ref. [29] . 

Far away from a critical point, Eq. Q is valid, so C 
and Ps can only be functions of the dimensionless ratio 
B2I Jvi- However, it can also be shown that B2/J12 — 
J2 C /2Ji2 = 7/2 for the B c = transition [3SJ, leading to 

C = 1/^1 + 7 2 /4, Ps = (l + C)/2, (7) 

in this limit. We have also solved for the critical behavior 
in more general cases, as described in [28], although the 
expressions are more complicated. 

In the limit of large chain size, JV > 1, we have pre- 
viously shown that 7 — > as ~ A -1 / 2 [5] [31] • Hence, 
C and Ps approach the constant value of 1, and the sin- 
gularity is suppressed, as consistent with Fig. [3J Using 
Eqs. (|6| and Q, we can define the depths of the C and 
Ps valleys as the difference between their asymptotic val- 
ues at large and small B c . In the limit of 7 — > 0, we find 
that the C and Ps valley depths approach zero quickly, 
as 7 2 /8 - 1/A. 

Triple-dot experiment. — We now propose an experi- 
ment to investigate the QPT in the opposite limit, N = 1. 
We consider the triple quantum dot geometry shown in 
the inset of Fig.[5ja). For simplicity, we assume the dots 
are singly occupied. The double quantum dot (on the 
left) is used to probe a spin "chain" of length 1 (on the 
right), whose ground state properties change dramati- 
cally as a function of the applied field B c at the transition 
point B c = 0. The field is applied only to spin c, neces- 
sitating a gradient scheme to cancel out the field on the 
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probe dots. The experiment proceeds by first preparing 
the triple dot in its ground state. We then turn off the 
exchange coupling to dot c and detune the probe double 
dot, so that the electron in dot 2 moves to dot 1 only 
when the probe is in a singlet state, due to the large 
singlet-triplet energy splitting [21]. The singlet proba- 
bility can then be monitored using charge sensing, via a 
nearby charge sensor. 

For the case N = 1, it can be shown that £>2 = ± J2 C /2 
m Eq. @ [25], so 7 = J 2c / 'Ji2- We can then determine 
the depth of the C and Ps valleys, as well as their widths, 
as a function of this ratio, as shown in Fig. [5] Generally, 
we see that applied fields of order B c ~ J12 are needed, to 
observe the QPT. Since quantum dot exchange couplings 
are typically of order /xeV, the magnetic field differences 
between the quantum dots needed in this experiment are 
of order 10 mT, and relatively easy to achieve in the 
laboratory 32J. 

Summary and conclusions. — In summary, we have 
shown that QPTs in a quantum dot spin chain cause non- 
analytic behavior in the entanglement between qubits at- 
tached to the chain locally or nonlocally, and we have ar- 
gued that there are benefits to using a local probe. Based 
on these methods, we proposed an experimental protocol 
for observing a QPT in a triple quantum dot. Although 
we do not describe the results here, we have also found 
that QPTs can be observed in other types of spin chains 
{e.g., XY), using a local probe. We believe that the local 
probe could be very useful for investigating entanglement 
in many quantum dot geometries, because of the practi- 
cality of performing singlet projection measurements. 

This work was supported by the DARPA QuEST pro- 
gram through a grant from AFOSR, and by NSA/LPS 
through grants from ARO (W911NF-08- 1-0482 and 
W911NF-09-1-0393). 
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SUPPLEMENTAL MATERIAL 



Here we provide detailed derivations of analytical expressions for the concurrence C and the singlet probability Pg 
near the critical transition at B c = in an odd-size chain, as reported in the main text. For completeness, we also 
provide expressions for the effective couplings J and local fields B, as originally derived in Ref. [T7] 



Expressions for effective couplings and local fields 



The following formulae, used in the main text, were first derived in Ref. |17j . 
For the nonlocal probe in the noncritical regime [Eq. (1) of the main text]: 



J\1a 



m>0 



J\cJl 



B x = Ji c (0|si z |0) , 

B 2 = J2c(0| SjVz |0) , 

\s la \m)(m\s Na \0) (for a = x, y, z) 



£o 



(1) 
(2) 

(3) 



Here \m) corresponds to the m-th eigenstate of the spin chain Hamiltonian, H c , with energy e m . Here, m — 
corresponds to the unique ground state. 

For the nonlocal probe in the critical regime [Eq. (2) of the main text], we define the pseudospin "up" state 
and "down" state (l-IJ-)) a s the ground states of the spin chain on either side of the transition, for higher and lower 
magnetic fields, respectively. The effective variables are 



B-2 -- 



2 

~2~ 



\siM) , 

\snM)) , 



J\2a — —J\cJlc 

J lea. 



B c — £4 - e-ff , 
SiQ,|m)(m|sjv a |-f|-) (V\\s\o\m)(m\s Na \\) 

^ic(tr|si+l^) (fortt = i,j) 

Jic((H\ai,\H) - (M\8 U m ({ora = z) , 

J2c{1\\sn+W) (fora = x,?/) 

J 2c «-fr|sjv a |lr> - (WI«J\r»W» (for a = z) . 



(for a = x,y,z) , 



For the local probe in the noncritical regime [Eq. (4) of the main text] 

B 2 = J 2c (0|si z |0) . 
For the local probe in the critical regime [Eq. (5) of the main text] 

•he 



B 2 = --f({t\sum 



J2ca — 



B c — £4 - % , 
J2c(1t\si+\ty) ({ora = x,y) 

J*c({1t\'u\1t) ~ m*u\H)) (iova = z) . 



(4) 

(5) 
(6) 
(7) 

(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



Concurrence of a nonlocal probe 

In this section, we derive Eq. (3) of the main text, which describes the concurrence of a nonlocal probe coupled to 
a pseudospin in the critical regime near the B c — transition. The pseudospin effective Hamiltonian is given by 



^cp — «/(Si + S2) ■ S c — B C S C:Z 



(14) 
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Without loss of generality, we will only consider B c > 0. S to t,z = S\ >z + S2,z + S c . z and Sf 2 = (Si + S2) 2 are good 
quantum numbers. For the subspace with Stot,* = ±3/2 and £12 = 1, the eigenstates are \Si iZ , S c . z , S2, z ) = | ttt) 
and I Hi), respectively. The eigenenergies are £±3/2 = J/2 T B c /2. For S to t,z = 1/2 and 6*12 = 0, the eigenstate is 
(I ttl) — I itt))/V2 and the eigenenergy is ~B c /2. For S to t, z = 1/2 and S12 = 1, the subspace is two-dimensional, 
and in the ordered basis {| tit), (I ttl) + I ltt))/V2} we have 



= -J 



J 



/2 



v^2 



2 J 2 



(15) 



The lowest eignevalue in this subspace is £1/2 = — J/4 — Ai/ 2 /2, where A1/2 = y § J 2 + — 

and S12 = 0, the eigenstate is (| til) — I IIDVv^ and the eigenenergy is B c /2. For 5 to t,z = 
subspace is two-dimensional, and in the ordered basis {(| til) + | lltiO/v^, | |tl), } we have 



JB C . For Stat,* = -1/2 
-1/2 and S 12 = 1, the 




(16) 



The lowest eignevalue in this subspace is E_ 1 / 2 = — J/4 — A_ 1 / 2 /2 1 where A_ x / 2 = y | J 2 + £? 2 + J^c- 
When B c > 0, we can easily see that E_ 1 / 2 is the ground state energy, and its eigenstate is given by 



I*gs) = -= ^cos - 1 til) + sin -| ||t) 



(17) 



where 9 is defined in the range < 9 < tt by 



cos6> 



sin6> 



1 + B c 
-4-1/2 

V2J 



/2 



(18) 
(19) 



We can construct the eight-dimensional density matrix for the ground state and trace out the pseudospin degree of 
freedom. The resulting reduced density matrix for the probe qubits is then given by 



P : 



icos 2 ^(|tl)(tl 



+ I It) Ut I + I tlXlt I + I ItXtl I) + sin 2 -I U)(U I 











I cos 2 f \ cos 2 § 
I cos 2 § \ cos 2 § 



\ 

















(20) 
(21) 



The concurrence is defined as [23] C = max{0, */ei — ^/e2 — \/^3 — -y/^t}; where E\ > e 2 > £3 > £4 > are the 
eigenvalues of the operator R = p (a y <g> a y ) p* (a y ® a y \ and a y is a Pauli spin matrix. From Eq. ( 21 ) , we then obtain 



C = cos 2 - = 



1 



1- 2 



B r 



1/2 



1 - 



B r 



\.P + Bl + JB C 



(22) 



Concurrence and singlet probability of a local probe connected to an odd-size spin chain 

In this section, we first derive Eq. (7) of the main text, which describes the concurrence and singlet probability of 
a local probe coupled to a pseudospin in the noncritical regime near the B c = transition. From Eq. (4) of the main 
text, 

H nc — B 2 S 2yZ + J12S1 • S2 . (23) 
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<Stot,z is a good quantum number. We readily obtain the two eigenstates | tt)i I 4-4-)) corresponding to Stot,z = ±1 
with eigenvalues E ±1 = J i2 /4± B 2 /2. In the subspace of S to t,z = 0, the Hamiltonian matrix in the basis {| t4-)> I It)} 
is given by 



Hn 




Jl2 

2 



-J 12 



(24) 



The lowest eigenvalue in this subspace is Eq = — Ji 2 /4 — \jJ\2^B\j2 and the eigenstate is l^o) = cos|| ^4- 
) — sin 1 1 Y\) 1 where 9 is defined in the range < 9 < tt by 



cos 9 



1 



1 + B; 



sm» = 



1-B 2 /JJ: 



J 2 
'12 



(25) 



(26) 



We can easily verify that |*o) is the ground state for any value of B 2 . From this, we obtain the expressions for the 
concurrence and singlet probability given in Eq. (7) of the main text. 

Next, we derive general equations for the concurrence and singlet probability of a local probe coupled to a pscudospin 
in the critical regime of the B c = transition. Equation (6) is a special case of this result. 

The pseudospin effective Hamiltonian in this case is given by 



Hep = J12S1 ■ S2 + J2cS 2 • S c 



B c S c ,z 



(27) 



Stot.z = Si tZ +S2.z + S c . z is again a good quantum number and we divide the full Hilbert space into subspaces with fixed 
Stot.z values. For Stot,z = ±3/2, there are two one-dimensional subspaces with eigenstates \Si^ z , S 2 . Zl S c _ z ) = \ ttt) 
and I Hi). The eigenenergies are £±3/2 = (J\ 2 + </2c)/4 T B c /2, respectively. For Stot,z = 1/2, the subspace is three 
dimensional and the Hamiltonian matrix in the ordered basis { | ttl) ; I tit) ; I ttt) } is given by 



H = 



/ J 12 -.hc+2B c 
1 A 

2 





2 

- J12 — Jlc — 2B C 
4 
J12 
2 





J12 
2 

- J12 + ^cz 2B ( 
4 



(28) 



We can obtain the eigenvalues by solving the cubic function \H — XI\ = 0. The ground state energy in this subspace 
is 



Ai = I J12 + J 2c + 2B C + 2V-D1 cos J 



where <j) is defined in the range < < 7r by 



Do 



COS ( 



sm(f> - 



2D-T 
V27A 



2D 



3/2 



and 



(4Ji 2 5 



4J. 



2c 



ABi -J, 2 J: 



12-J2c 



4Ji2-B c ^ 2J2 C B C 



D 2 = 16 (8 



8J 12 + 8J 2c — 8B C — 3J 12 J2C — 3Ji2 J 2c + 6J12 J2cB c — \2J\ 2 B C + 6J2 C B C + \2J 12 B C — 6J 2c B, 



A=^(ADl-Dl) . 



(29) 

(30) 
(31) 

(32) 
(33) 
(34) 



Without loss of generality, we will only consider B c > 0. In this case B c > 0, the ground state energy Ai in the 
Stot,z = 1/2 subspace is always lower than the ground state in the subspace with S to t,z — ±3/2. The corresponding 
eigenstate of the three-spin system is 



|*GS>=*l|m)+Z2|Ut> + Z3|ltt) 



(35) 
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where 



zi = - 



Zl = 



z 3 



Z = 



2y/Z(F-X 1 ) ' 
1 

7z ' 

Jl2 

2y/Z(F + X 1 ) ' 

(%) 2 (Ax + Ff + (^) 2 (Ai - Ff + (A? - F 2 )' 



(\i-F>y 



Jl2 ~ Jlc + 2i? c 

4 



(36) 
(37) 
(38) 

(39) 

(40) 



As before, we construct the eight-dimensional density matrix for the ground state and trace out the pseudospin degree 
of freedom. The resulting reduced density matrix for the probe double dot is then given by 



p = z\\ ttxtt i + z\\ nxn i + 4\ n>ut i + ^ (i n)ut i + 1 itxn d 

/ zf \ 

z| z 2 Z3 

z 2 z 3 z 2 

\ / 

We can then compute the concurrence, which is given by 

C = 2\z 2 z 3 \ 



(41) 
(42) 



J12 



Z\F + \A 



(43) 



The singlet probability Ps is given by 



Ps = (%|5) = 2^2-^ 3 ) 2 



1 

2Z 



1 - 



J12 



2(F + Ai) 



(44) 



We could obtain Eq. (6) of the main text by plugging B c = into these expressions. However, the algebra is rather 
complicated. Here, we solve the effective Hamiltonian specifically for the case B c = 0: 



H — J12S1 • S 2 + J2c§2 • S c . 

Similar to the nonlocal probe case, we obtain the ground state 



1 6 1 . 6 

—= cos - H ■= sin - 

V2 2 ^6 2 



|*G5> : 

where 6* is defined in the range < 9 < tt by 

cos# 



tit) - I ltt>] 



3 sin -I m) 



J11 



J2c 



sin# = 



V3 7- 
-9- <J2c 



J 12 — J12J2C + >^2c 



(45) 
(46) 

(47) 
(48) 



Calculating the reduced density matrix after tracing out the pseudospin degree of freedom, we obtain Eq. (6) of the 
main text. 

I 



